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Abstract. Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad- 
band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious 
similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from 
individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact 
feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this 
hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars 
using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling. 



1. Introduction 

Blazar jets are known to radiate across the entire elec- 
tromagnetic spectrum (e.g. Marschcr 2006, I2010p . Their 
radio to UV (sometimes - X-ray) emission is believed 
to be dominated by synchrotron radiation of relativis- 
tic electrons while radiation at higher energies could be 
due to the inverse Compton scattering of synchrotron 
photons emitted by the electrons themselves (the syn- 
chrotron self-Compton process, SSC, IJones et al. 19741 
Ghiscllini & Maraschi 1989|) and photons from exter- 
nal sources (External Compton, EC, Sik ora et al. 1 994, 
IDermer & : Schlickcisc r 2002[ ). The sources of the external 
seed photons for the EC process include the accretion disc, 
broad line region (BLR) clouds, warm dust and the cosmic 
microwave background (CMB), with their relative contri- 
bution vary for different blazars. 

Despite the wide acceptance of the picture outlined 
above, the exact location, geometry and physical proper- 
ties of the regions responsible for blazar emission at dif- 
ferent bands remain largely unknown. Considerable suc- 
cess was achieved by single-zone models (assuming a sin- 
gle spherical relativistically moving blob of magnetized 
plasma which emits synchrotron radiation and interacts 
with an ambient photon field) in explaining blazar emis- 



sion from IR to 7-ray energies (among the many recent ex- 
amples: lAbdo et al. 20101 IFinke fc Dermer 20101 and the 
references in Table [TJ. A significant fraction of blazar 
radio emission is known to originate at the extended 
parsec-scale jet and cannot be accounted for in the frame- 
work of a single-zone model. However, emission from in- 
dividual jet components resolved with the Very Long 
Baseline Interferometry (VLBI) can usually be well de- 
scribed by a single uniform synchrotron source model 
(Pacholczyk 1970). It is tempting to identify one of the 
jet components resolved by VLBI with the zone responsi- 
ble for the emission at higher frequencies. 

In this contribution we describe a program of coordi- 
nated multi-frequency (4.6-43.2 GHz) Very Long Baseline 
Array (VLBA), Swift and Fermi observations of selected 
7-ray bright blazars with the aim to identify the possible 
source of high-energy radiation within the parsec-scale jet 
and put constraints on physical properties based on radio 
spectra. 

2. Observations and analysis 
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We have used the VLBA (Napier 1994) to observe twenty 
blazars that were expected to be bright 7-ray emit- 
ters, prior to the launch of Fermi-GST. The observations 
were conducted simultaneously at seven frequencies (4.6- 
43.2 GHz) utilizing the array's unique capability of fast 
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Table 1: Comparison of the emission region parameters 
estimated from the VLBA data and SED modeling 



Multi-frequency VLBA 


1 Published SED model 


results (this work) 


parameters 



R43 GHz PvLBA ^VLBA. D va , T e f SED -RsED PsED &SED Ref. 

[10 1 " 



[G] 



[10 11 



AO 0235+16 


< 5800 


0.8 


< 11 


24.0 












B0528+134 


< 7600 


1.4 


< 1.2 


31.2 












S5 0716+714 


< 2100 


a = 


0.4 d 


10.9 


11 


40 


2.0 


1 


(1) 


OJ 248 


< 6300 


1.6 


< 23 


13.1 












OJ 287 


< 2800 


a = 


0.7 d 


17.0 












W Com 


< 1700 


0.8 


< 118 


1.2 


20 


3 


2.55 


0.35 


(2) 


3C 273 


< 1800 


2.0 


< 0.2 


17.0 


9 


20 


2 


12 


(3) 


3C 279 


< 4600 


1.4 


< 14 


24.0 


21.5 


25 


2.0 


1.8 


(4) 


PKS B1510-089 


< 3900 


a = 


0.2 d 


16.7 


37 


18 


1.9 


0.09 


(5) 


4C 38.41 


< 3500 


1.0 


< 1.8 


21.5 












Mrk 501 


< 470 


1.6 


< 67 


f 


20 


1.03 


2 


0.3 


(6) 


NRAO 530 


< 6200 


1.6 


< 20 


10.7 












B1959+650 


< 510 


1.4 


< 300 


f 


18 


7.3 


2 


0.25 


(7) 


B2155-304 


< 2400 


1.0 


< 745 


f 


32 


150 


1.3 


0.018 


(8) 


BL Lac 


< 590 


1.0 


< 3 


7.3 












3C 454.3 


< 6500 


a 


0.8 d 


33.2 


26 


15 


1.1 


5.4 


(9) 


B2344+514 


< 680 




< 57 c 














a p is the power 


law index 


■ the electron 


energy 


distribut 


ion N(E) = 


NoE-" 


For 



the optically thin part of the synchrotron spectrum p — 1 — 2a where a is defined as 
Si, ~ v a . The estimates correspond to the parsec-scale radio core. c The value is in 
the observer's frame. The homogeneous synchrotron source model is not applicable 
for such spectrum. e The variability Doppler factor from Hovatta et al. (2009) used 
to transform the observed magnetic held strength into the source frame. The same 
Doppler factor as in the corresponding SED model was adopted. 
Refe rences: l:|Chen et al. (2008)| 2: |Acciari et al. (2009)| 3: |Pacciani et al. (2009)1 

'Aiiderhub et al. (2009*7] 7- 



4: |Giuliani et al. (2009)| 5: ~ this work; 6 : 
ITagliaferri et al. (2008)| 8: |Aharonian et al. (2009)| 9: 



Bonnoli et al. 



(2010)| 



frequency switching between individual VLBI scans. After 
the initial calibration in AIPS (Grciscn 1990), the sources 
were self-calibrated and imaged independently at each 
frequency using Difmap (She pherd 1994[ ). This software 
was also used to model the source structure in the uv- 
plane and derive constraints on the core size at 43.2 GHz 
(Tabled]). 

A special procedure was developed to improve am- 
plitude calibration of the correlated flux density result- 
ing in a ~ 5 % accuracy in 4.6-15.4 GHz range and a 
~ 10 % accuracy at 23.8 and 43.2 GHz. Images at dif- 
ferent frequencies were aligned with each other using op- 
tically thin regions of the jet to compensate for phase 
center shifts and possible frequency-dependent core po- 
sition ("core shift", ILobanov 19981 IKovalev et al. 2008j) . 
We applied a spectrum extraction technique described by 



Sokolovsky et al. (2010) which relies on the positions and 
flux densities of individual CLEAN components as op- 
posed to the use of restored VLBI images in order to sup- 
press spectral artifacts arising from the convolution with 
a Gaussian beam. 

The sources were observed with Swift in X-ray and 
optical-UV bands typically within two days after the 
VLBA observation. The Swift data were analyzed us- 
ing standard tools provided in the HEAsof^ package. 
Whenever feasible, standard stars within the UVOT field 
of view were utilized to improve photometry in UBV 
bands. The Galactic correction was applied using the val- 
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Fig. 1: Quasi-simultaneous SED of PKS B1510-089 
constructed using observations with the VLBA, 
Swift, Fermi/LAT, NOT, SMA and the 2.1 m tele- 
scope Guillcrmo Haro. The violet curve (a) represents the 
combined SED model: green curve (b) is the synchrotron 
component, red (c) is the SSC component, blue (d) is the 
accretion disk and brown (e) represents EC scattering of 
the disk radiation (see also Table [2]) . 



ues from |Schlegel et al. (1998)| and |Kalberla et al. (2005) 
for the optical-UV and X-ray absorption respectively. 

To obtain the quasi-simultaneous SED of 
PKS B1510-089 (Fig. P discussed below, the VLBA, 
Swift and Feum/LAT data collected during two days 
2009 April 09 - 10 were complemented by UBVR c i 
photometry obtained on April 12 with the 2.56 m Nordic 
Optical Telescope (NOT) and infrared JHK S photometry 
obtained on April 07 and 17 with the 2.1 m telescope 
of Guillermo Haro Observatory in Cananea, Mexico. 
The Submillimeter Array (SMA) provided flux density 
measurements at 1 mm from 2009 April 05 and 14. 
These observations resulted in the most well-sampled 
quasi-simultaneous SED (containing multi-frequency 
VLBI data) among all our sources, and for that reason it 
was the first we chose to construct a model of. 

UBV magnitudes of PKS B1510-089 observed by the 
Swift/UVOT are systematically (0™5-0™3) brighter then 
those observed by the NOT two days later. Since both 
observations were calibrated against the same set of com- 



1 http : //heasarc .nasa.gov/lheasof t/ 



parison stars from Gonzalez-Perez et al. (2001) the differ- 
ence can be interpreted as a clear sign of optical inter-day 
variability. 

3. Discussion of the first results 

Parsec-scale radio emission of all the observed sources 
(Table[T]) is dominated by a bright unresolved core (the ap- 
parent origin of the jet, see a discussion by Marschc r 2008[) 
which exhibits high-amplitude flux density variability (e.g. 
IVercellone et al. 2010|) . The parsec-scale core is a natu- 
ral candidate to be directly related to the bright variable 
emission at 7-rays (Kovalcv et al. 200"9]) and other bands. 
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Power law fit: 
a = 0.2 



Power law fit: 
a = 0.0 



Synchrotron fit; 

S peak = 1 ' 24 ^ 

a = -0.3 




Frequency (GHz) 

(a) PKS B1510-089 



Frequency (GHz) 

(b) OJ 248 



Fig. 2: VLBA spectra of core regions of PKS B1510-089 
(a) and O J 248 (b) . The former spectrum is fiat while the 
latter one shows a prominent synchrotron self- absorption 
peak. The green line is a power law fit. The blue curve is 
a homogeneous synchrotron source model. 



The multi-frequency VLBA observations have detected 
the synchrotron self-absorption peak in core regions of 
most of the observed sources with the exception of 
B0716+714, OJ 287, PKS B1510-089 and 3C 454.3. The 
latter sources show inverted spectra in the 4.6-43.2 GHz 
range, characteristic of a partially optically-thick inhomo- 
geneous source. These sources may have a detectable self- 
absorption turnover at frequencies > 43 GHz. Fig. [3] com- 
pares core radio spectra of PKS B1510-089 and OJ 248. 

For sources where the self-absorption peak was de- 
tected, it was possible to put upper limits on the mag- 
netic field strength using the procedure described by 



Sokolovsky et al. (2010) The source intrinsic values were 
computed using variability Doppler factor estimations 
from |Hovatta et al. (2009) The resulting liberal upper 
limits on the magnetic field strength owe to poor reso- 
lution at lower observing frequency (4.6 GHz) which ulti- 
mately limits our ability to unambiguously extract spec- 
trum of a given spatial region. 

In the case of a bright isolated component, which dom- 
inates the emission across the whole observed frequency 
range, it is possible to estimate its size using the high- 
est observing frequency (43.2 GHz) instead of the lowest 
one and, therefore, obtain much tighter constrains on the 
component size and the magnetic field strength. The ap- 
plicability of this "isolated component" scenario to each 
individual source is currently being investigated. 

Table [T] presents constrains on size (R43 gh z ), mag- 
netic field strength (-Bvlba) and electron energy spec- 
tral slope (pvlba) in the radio core region, see 



Sokolovsky et al. (2010) for a more detailed discussion. 
These values are compared to the corresponding parame- 
ters of SED models (i?sED, -Bsed, Psed) reported in the 
literature. Doppler factors from SED models (L>sed) and 
those assumed for -Bvlba computation (D mr ) are also 
listed. The SED model parameters are well within the con- 
straints set by the multi-frequency VLBA data, except for 
3C 273 which shows a discrepancy in -Bvlba _ -Bsed that 
may be due to insufficient resolution of our data set. There 



is evidence from 86 GHz VLBI data that B ~ 1 G in the 
mm-wavelength core of 3C 273 ( Savolainen et al. 2008). 

To further test the possible relation between the radio 
core and the emitting blob implied by single-zone SED 
models, we have constructed the quasi-simultaneous SED 
of PKS B1510-089 presented in Fig. □ The use of the 
VLBI core flux density (instead of the total flux density 
obtained by single-dish radio observations) allows us to 
exclude the emission coming from the parsec-scale jet and 
larger scales. Our ambition was to check how this will 
affect the SED modeling. 

Table 2: PKS B1510-089 SED model parameters 



Minimum Lorentz factor 


7niin — 


= 30 




Brake Lorentz factor 


7brk = 


= 1.9 x 10 3 




Maximum e~ Lorentz factor 


7max — 


= 1.0 x 10 5 




e~ energy slope below 7 br k 


Pi = 


= 1.9 




e~ energy slope above 7 br k 


£>2 = 


= 3.9 




Doppler factor 


Z> = 


= 37 




Bulk Lorentz factor 


r = 


= 37 




Magnetic field strength 


B = 


= 0.09 G 




Variability time 


^var — 


= 2.16 x 10 4 


sec 


Blob radius 


i? = 


= 1.8 x 10 16 


cm 


Jet power (magnetic field) 


Pjet, B = 


= 2.8 x 10 43 


erg/s 


Jet power (electrons) 


P. 

jet e 


= 7.5 x 10 45 


erg/s 


Black hole mass 


M BH = 


= 1 x 10 9 M Q 


Accretion efficiency 


77 = 


= 1/12 




Eddington ratio 


idisk/LEdd = 


= 0.06 




Blob distance from disk 


T = 


= 5.1 x 10 17 


cm 



Contribution from a hot thermal component, proba- 
bly the accretion disk, is evident in the optical and UV 
bands. The SED in the IR to 7-ray range can be ex- 
plained by the single zone EC model with the Compton 
scattered accretion disk photons being responsible for the 
bulk of the observed 7-ray emission. In this model, the 
blob is close enough to the accretion disk so the accre- 
tion disk photons can reach the blob directly, with no 
need to be scattered back from BLR clouds or dust. Note, 
that due to rclativistic aberration some of the disk pho- 
tons will be coming to the blob (nearly) head-on. The 
photons will primarily come from a radius of the disk 
of y/3r where r is the blob's distance from the black 
hole (|Dermer fc S chlickcis er 2002| . The blob size was con- 
strained by the observed 7-ray variability time scale of the 
order of a few hours (jAbdo et al. 2010)) . Parameters of the 
model are summarized in Table [51 We note that the model 
fit is not unique, single-zone models with a different source 
of seed photons could probably also produce an acceptable 
fit (cf. lAbdo et al. 20101) . 

The constructed model (Fig. [T]) represents well the IR 
to 7-ray data, however it dramatically underpredicts the 
radio emission because of the synchrotron self- absorption 
occurring in the small blob at a high frequency. A sig- 
nificantly larger blob size would be inconsistent with 
the observed variability timescale. Even if the variabil- 
ity timescale argument were dismissed, a single-zone SED 
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model with a larger blob would have difficulties in ex- 
plaining the observed hard radio spectrum which is incon- 
sistent with the softer X-ray spectrum or optically thick 
synchrotron radiation from a uniform source. One way to 
overcome this difficulty would be to introduce a second 
brake in the electron energy spectrum. 

Alternatively, one could abandon the attempt to ex- 
plain the whole SED with a single-zone model and assume 
that even the radio emission observed from the cm-mm 
band core is coming from a larger (-Rsed < R < R43 GHz) 
structure downstream of the blob responsible for the IR- 
to-7-ray emission. This structure (which cannot be ob- 
served separately from the blob due to limited resolu- 
tion of the available VLBI data) could be a smooth jet 
(Blandford & Konigl 1979) or a number of distinct jet 
components (|Marscher 1980[) - perhaps the blobs which 
contributed to high-energy emission earlier... One could 
argue that currently this extended region does not con- 
tribute significantly to the optical and 7-ray emission on 
the basis of the observed short timescale variability in 
these bands. A more detailed SED modeling is needed to 
test this scenario. 

4. Summary 

The program of coordinated multi-frequency (4.6- 
43.2 GHz) VLBA, Swift and Fermi observations of selected 
7-ray bright blazars is described and the first analysis re- 
sults are discussed. Constraints on the size (R43 GHz) and 
magnetic field strength (Bvlba) in the parsec-scale core 
region derived from the VLBA observations are consistent 
with the values (Rsed, -Bsed) usually assumed in single- 
zone SED models (Table[T]). However, the single-zone SED 
model which we tested for PKS B1510-089 has difficul- 
ties in explaining the observed level of core radio emission 
while being consistent with the observed 7-ray and opti- 
cal variability timescales. The multi-frequency VLBA and 
SED data analysis for the whole sample is in progress. The 
SED data on PKS B1510— 089 discussed here are available 
from KVS upon request. 
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